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■ Abstract 



. We establish the applicability to transport phenomena in the cuprate super- 

conductors of a nearly antiferromagnetic Fermi liquid (NAFL) description of 
, the magnetic interaction between planar quasiparticles by using it to obtain 

Q\ , the doping and temperature dependent resistivity and Hall conductivity seen 



C3 



experimentally in the normal state. Following a perturbative calculation of 
the anisotropic (as one goes around the Fermi surface) quasiparticle lifetimes 



C. 

O , which are the hallmark of a NAFL, we obtain simple approximate expressions 

O 

for the longitudinal, o~ xx , and Hall, a xy , conductivities which reflect the mag- 
netic crossovers seen experimentally as one varies the doping level and temper- 
ature. We present a simple phenomenological model for the variation in mean 
free path around the Fermi surface, and use this to extract from experiments 
on a xx and a xy quasiparticle lifetimes in the hot (strongly coupled quasipar- 
ticle) and cold (weakly coupled quasiparticle) regions of the Fermi surface 
which are consistent with the perturbation theory estimates. We improve 
upon the latter by carrying out direct numerical (non-variational) solutions 
of the Boltzmann equation for representative members of the YBa2Cu306+z 
and La2_ x Sr x Cu04 systems, with results for transport properties in quan- 
titative agreement with experiment. Using the same numerical approach we 
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study the influence of CuO chains on the a-b plane anisotropy and find results 
in agreement with experimental findings in YE^CuziOs. 
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I. INTRODUCTION 



The peculiar normal state properties of High Temperature Superconductors (HTS) have 
been extensively studied over the past decade in the expectation that their understanding 
will likely reveal the mechanism of superconductivity. Yet no consensus has emerged as to 
their physical origin. Indeed, it is striking that so much work has been done on understanding 
a variety of complicated experimental situations, while rather basic measurements, such as 
the d.c. transport, are not well understood. So far it has been clearly established^ and well 
documented in the literatureH that practically all normal state properties in these materials 
are anomalous with respect to the behavior of standard Landau Fermi Liquids (FLs). This, 
in turn, poses the question of whether Landau FL theory is applicable to HTS.i A variety 
of proposals based on non-Fermi Liquid©"! (NFL) or marginal FlI theories have emerged 
which have had considerable success in explaining a number of specific experimental results. 
However, at present no single theory has been able to account for all the anomalies found 
in the normal state properties of the cuprates. 

In this paper we focus on the planar resistivity and the Hall effect, which have drawn par- 
ticular attention in the HTS community, due to their arguably "contradictory" behavior!^ 
The temperature dependences of the Hall angle, 6h, and the longitudinal resistivity, p xx , 
are so vastly different that the Hall coefficient, Rh, in practically all cuprates is a strong 
function of temperature, especially near optimal doping. This phenomenon by itself, is 
rarely observed in ordinary FLs.0 It is also found that all HTS near optimal doping have 
p xx proportional to the temperature, T, in a wide region of interest. Hence it appears that 
the transport properties in the presence of magnetic field as well as in its absence require 
two separate scattering mechanisms, and presumably two separate relaxation rates. In fact 
this is the only common scheme among the competing theories for the explanation of the 
anomalous transport in HTS materials. 

We begin with a brief review of the theoretical proposals. Anderson, who first drew 
attention to the problem of explaining the unusual magnetotransport in the normal state,0 
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has proposed a scenario based on a NFL ground state. In his scenario the two rates observed 
in experiment are attributed to scattering of two distinct many-body excitations, spinous 
and holons. The existence of these excitations is well established in ID Luttinger liquids 
and Anderson has put forth a non-trivial generalization of this concept to the quasi-two 
dimensional system found in cuprate superconductors. Strictly speaking one can talk about 
a NFL in any system at, or near half-filling, which has an interaction sufficiently strong 
to produce the Neel instability, but the validity of this approach at large doping levels is 
somewhat questionable. In the Anderson scenario both spinons and holons contribute to the 
resistivity, but it is often tacitly assumed that the scattering of holons dominates, while only 
spinons, which effectively commute with the magnetic field, B, appear in the temperature 
variation of the Hall angle. Although commonly used to analyze experiments, this proposal 
has yet to be verified by an explicit calculation of the transport coefficients.lll 

Coleman et alHH have developed a phenomenological transport theory based on the as- 
sumption that the dominant scattering in cuprates is mediated by an interaction which 
distinguishes quasiparticle wave functions with respect to their charge conjugation symme- 
try. Such interactions are absent in ordinary FLs and hence this can be regarded as a NFL 
model, although in this case there is no spin-charge separation and hence the model is more 
correctly categorized as an unusual FL. In this model quasiparticles with different symmetry 
have vastly different scattering rates, which in turn have different temperature dependences, 
ultimately leading to the correct behavior of all transport coefficients. Although the origin of 
the unusual interaction has not been identified, this is the only model to date which appears 
to be consistent with the recent experimental results on optical reflectivity.0 

The anomalous Hall effect can also be understood in terms of the bipolaron theory, devel- 
oped by Mott and Alexandrov.0 Anomalous behavior in their model arises from Anderson 
localization of bipolarons: at low temperatures bipolarons in a random potential lose their 
mobility, producing an enhanced value of Rh- At sufficiently high temperature, larger than 
the random potential binding energy, the mobility is recovered, and Rh saturates, as seen 
in experiment. The bipolaron model is also consistent with a number of photoemission 



experiments. 

A variety of FL-based models have also been utilized to explain the anomalous trans- 
port in cuprates. Clearly, in order to account for experiments, one must introduce the two 
scattering rates without violating the FL ground state. This is done either by introducing 
a certain anisotropy of scattering or by assuming a different mechanism of scattering in 
the presence of magnetic field. Initially, the anomalous Hall effect was attributed to the 
well-known skew scattering present in many materials with strong spin-orbit coupling]^ 
the anomalous temperature variation of p xy arises from anisotropic scattering of conduction 
electrons by local magnetic moments. This approach naturally produces the observed tem- 
perature dependence of p xy , provided the resistivity is assumed linear in T. However, one 
runs into problems for systems with doping levels for which the resistivity is not propor- 
tional to temperature. Moreover, skew scattering should saturate with increasing magnetic 
field, due to the complete polarization of localized moments. Such saturation has not been 
observed, although some discrepancy from the linear in B Hall resistivity does occur in 
extremely high pulsed fields.0 

Levin and QuadeiH focused on a very important subgroup, the bilayer cuprates. In 
their approach the anomalous behavior of the Hall effect is due to the multiband structure 
of bilayer cuprates. Their main assumption is that one of the bands is very close to the Fermi 
level and responsible for the anomalous resistivity, while the other is responsible for the Hall 
angle. The former is thermally activated at higher temperatures in underdoped materials 
so that the model can explain qualitatively the behavior of resistivity in the underdoped 
bilayers. Although their approach appears reasonable, the striking similarity between the 
experimental results on single layer and bilayer materials, which we discuss in Sec. [TV], 
indicates that this approach may have to be revised before it possesses general applicability. 

Another class of models is based on anisotropic scattering, which can arise from either 
anomalous band structure in the cuprates or spin-fluctuation scattering. Newns§ has pro- 
posed a model based on an anomalously large density of states near the Fermi surface (FS) 
for quasiparticles in the vicinity of the so called Van Hove singularities: angle resolved 



photoemission (ARPES) experiments in a number of samples show extended flat band fea- 
tures near k = (ir, 0) and symmetry related points in the Brillouin zone, which give rise 
to highly anomalous quasiparticle scattering. The two scattering rates necessary to explain 
the measured transport properties are those for quasiparticles near and far away from these 
singularities. Although the model can produce a temperature dependent Hall coefficient, 
Rh = p xy , it is unclear whether it can simultaneously account for the anomalous behavior 
of the resistivity. 

We next consider models based on spin-fluctuation scattering. One of the first con- 
tributions has been provided by Carrington, Cooper and their collaborators^ who used a 
phenomenological model of the spin-fluctuations and a realistic FS and calculated a num- 
ber of transport coefficients in cuprates. They assumed that large parts of the FS have an 
anomalous, linear in T scattering rate, while smaller parts have a Fermi liquid, quadratic 
in T scattering rate. Although the anomalous quasiparticle scattering rate, assumed to be 
due to spin-fluctuations, has been taken as an input to the theory, they obtained a remark- 
able agreement not only with the Hall effect, but with the thermoelectric power as well. 
This pioneering work motivated a number of authors to explore more particular choices of 
spin-fluctuation spectra and their relationship with transport 

Trugman@ calculated the band structure renormalization in the presence of strong an- 
tiferromagnetic (AF) fluctuations within the context of the t — t' — J model and obtained 
a single particle band structure which has the form of a precursor to a spin- density- wave 
(SDW) state. In this case one finds flat regions in the band, somewhat similar to those 
mentioned above, which are characterized by an energy scale which is much smaller than 
the fermionic bandwidth. In essence the size of the Fermi energy is effectively reduced to 
an energy fi<i, and hence, just as in ordinary FLs, the resistivity deviates from its usual 
T 2 dependence, while the Hall coefficient is temperature dependent. However we note that 
this anomalous temperature dependence of the Hall effect is primarily due to the strong 
anisotropy of scattering rates and densities of states in different regions of the Brillouin 
3. Similar models were also provided by Kim et al@ and by Dagotto et alS 



zone. 



Lercher and Wheatleyc3 also studied a spin-fluctuation model of magnetotransport in the 
cuprates. They concluded that the behavior of several quantities of interest is significantly 
altered depending on whether one incorporates Umklapp processes. At high temperatures 
they find a highly T dependent Hall coefficient, with a resistivity approximately linear in 
T. As in the above cases, their calculation is based on the anisotropy of scattering rates in 
different regions of the Brillouin zone, although here the anisotropy originates in the specific 
structure of the effective interaction. 

The nearly antiferromagnetic Fermi liquid (NAFL) model for transport in the cuprates 
resembles the model studied by Lercher and Wheatley, as well as the earlier Hubbard model 
calculations of the resistivity by Bulut, Scalapino and White,0 with one important differ- 
ence: the choice of the effective interaction. Although in both cases the effective quasiparticle 
interaction, V e ff, assumed to be of magnetic origin, is a strong function of momentum trans- 
fer, and is highly peaked near the antiferromagnetic wavevector Q, Lercher and Wheatley 
calculate V e ff using an RPA formalism with a local restoring force, while in the NAFL model 
one takes the momentum and frequency dependence of the interaction to be that of the spin- 
spin response function which provides a quantitative account of the NMR experiment s,il as 
discussed in the next Section. While an RPA formalism can provide insight into the physical 
origin of the spin-fluctuations, it neglects changes in the effective interaction brought about 
by the spin-fermion vertex corrections, which lead to a sequence of crossovers in the low 
frequency magnetic behavior in the normal state.il Thus it is not capable of accounting for 
the rich morphology of the magnetic properties seen in NMR experiments]! a morphology 
which we shall see is reflected in transport experiments. 

Pioneering calculations of the resistivity of YBa2Cu3C>7 using a NAFL model were car- 
ried out by Monthoux and Pines,0 who found, in a strong coupling (Eliashberg) calculation 
which incorporated vertex corrections, that a good quantitative account of experiment could 
be obtained using a spin-fluctuation spectrum taken from NMR experiments and a coupling 
constant which yielded a superconducting transition at 90K. This work provided an im- 
portant bridge between the anomalous normal state properties and the superconducting 
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transition. There remained the important question of whether an NAFL model could also 
account for the anomalous Hall effect. The present authors addressed this question by using 
standard Boltzmann transport theory to calculate both the resistivity, p xx , and the Hall 
conductivity, a xy of YBa2Cu3O7.il A direct (nonvariational) numerical solution of the BE 
showed that the highly anisotropic quasiparticle scattering at different regions of the Fermi 
surface (brought about by the highly anisotropic magnetic quasiparticle interaction which 
characterizes a NAFL) could give rise qualitatively, and in some cases quantitatively, to the 
measured anomalous temperature dependence of the resistivity and Hall coefficient while 
yielding the quadratic temperature dependence of the Hall angle. ME The present paper 
extends these calculations in a number of ways. First, we use perturbation theory to cal- 
culate, for a broad range of doping and temperature regimes, the anisotropic quasiparticle 
scattering rates around the Fermi surface; we then use these results to calculate the longi- 
tudinal and Hall conductivities. Our perturbation theory approach not only makes explicit 
the difference in quasiparticle behavior at hot spots (those regions of the FS for which the 
magnetic interaction is anomalously large) and in cold regions (the remaining parts of the 
FS), but enables one to examine in detail the changes in quasiparticle lifetimes in both 
regions brought about by both the crossovers in the spin-fluctuation spectra measured in 
NMR experiments and the changes in FS shape measured in experiments using angle re- 
solved photoemission techniques. On taking the relevant averages of quasiparticle lifetimes 
for the longitudinal and Hall conductivities, we are then able to obtain simple approximate 
expressions for these quantities as a function of doping and temperature which are in good 
agreement with experiment. 

The substantial anisotropy of quasiparticle mean free paths found in the perturbation 
theory approach led us to develop a phenomenological model for the variation in mean free 
path as one goes around the FS, and to use this to extract from experiments on a xx and a xy 
the behavior of the quasiparticle lifetimes, Thot and r co id, as one varies temperature and dop- 
ing level for the YBa 2 Cu30 6+;c (123), La2_ x Sr x Cu0 4 (214), and Tl-doped systems. Again, 
good agreement is found between the perturbation theory estimates and the experimental 
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results. 

Guided in part by the perturbation theory results and our phenomenological model, we 
then carry out direct numerical solutions of the BE for representative members of the 123 
and 214 systems, and for a highly simplified model of a system, such as YE^Ch^Os (124), 
in which the presence of chains brings about considerable planar anisotropy. For all three 
cases considered, good agreement is found between theory and experiment. 

The paper is organized as follows: in the next section we provide the necessary back- 
ground on the signatures of NAFL model in NMR experiments and define the relevant 



crossover temperatures. In Sec. Ill we analyze the properties of NAFLs using a simple 



perturbative approach. In Sec. [Iy| we review experimental results found for the 123, 214, 
and overdoped Tl systems, and call attention to the regularities (or lack thereof) which 
p xx (T) and p xy (T) exhibit as a function of doping and temperature. We then present our 
phenomenological analysis of a xx and a xy in terms of a highly anisotropic mean free path, 
and use this analysis to extract from experiment the quasiparticle lifetimes in the hot and 
cold regimes. In Sec. [V], we present numerical calculations for a realistic set of band and 
spin-fluctuation parameters, and compare our results to experiment. Finally, in Sec. [VT| we 
summarize our conclusions. 



II. BACKGROUND ON PHYSICAL PROPERTIES OF NAFLS 

In the NAFL description of the normal state properties of the superconducting cuprates 
it is the magnetic interaction between planar quasiparticles which is responsible for the 
anomalous spin and charge behavior. The magnetic properties of the system are specified 
by a dynamical spin-spin response function of fermionic origin, %(q, u), which near a peak 
at Qj, in the vicinity of the commensurate wavevector, Q = (tt, tt), takes the form: 

i 1 + (q - Qi)"V - iuj/u)sf - /A 2 
Here XQi ^ Xo is the magnitude of the static spin susceptibility at a (possibly) incommen- 
surate wavevector, Qj, near Q, £ is the antiferromagnetic correlation length, oj s f specifies 
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the low frequency relaxational mode, brought about by the near approach to antiferromag- 
netism, and A allows for the possible presence of a spin gap at large energies. We have 
assumed a system of units in which the lattice spacing a — 1. The quasiparticle spectrum 
is assumed to take the tight-binding form, 



where t, t' and t" are the appropriate hopping integrals, while the effective magnetic inter- 
action between the planar quasiparticles is specified by 



For a given system the parameters \Qn 6j ^sf which determine x(q, a>) are taken from fits 
to NMR and INS (inelastic neutron scattering) experiments, while the effective coupling 
constant, g, is assumed to be momentum independent for the wavevector near Q^. As 
discussed by Chubukov et al,S in principle the effective interaction, Eq. can be derived 
microscopically starting with, e.g., a one-band Hubbard model, but for a number of reasons 
it has not yet been possible to carry through such a derivation. 

Since in practice the effective damping of the magnetic excitations is large, the term 
proportional to u 2 in the denominator of Eq. (fj) can be neglected, so that we may write: 



efc = — 2t(cos k x + cos k y ) — At' cos k x cos k y — 2t" (cos(2k x ) + cos(2k y )) 



(2) 




(3) 




(4) 



UJ q — IUJ 



where we have introduced the scale factor, a, which relates \Qi to £ 2 , 



(5a) 




(5b) 



and 



^0 — ^s/£ ) 



(5c) 
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In the next Section we shall for the most part assume that the spin fluctuation spectrum 
possesses only a single peak at Q; we consider, in Sec. [V], the effects of incommensuration 
on the transport properties. 

The effective interaction V e ff, Eq. (|3]), has an obvious property: for sufficiently large 
correlation lengths it is highly peaked for momentum transfers in the vicinity of the anti- 
ferromagnetic wavevector Q. The importance of this fact is hard to overestimate: if the 
FS of the system of fermions, defined by the quasiparticle dispersion is such that it 
intersects the magnetic Brillouin zone Z (see Fig. |l|), then quasiparticles in the vicinity of 
these intersection points on the FS, often referred to as hot spot sP are much more strongly 
scattered by the spin- fluctuations than those which are on other parts of the FS. This is the 
main source of the anisotropy of quasiparticle behavior in this model, since the temperature 
variation of the quasiparticle scattering rates at and far away from hot spots is in general 
very different. We shall return to this point in the following Section. 

Barzykin and Pines (BP) have utilized Eq. ([I]) in their analysis^ of the NMR results in 
HTS. Because their work, subsequently extended to an analysis of INS results,il provides 
an important insight into the morphology of the magnetic spectrum in these materials, 
we briefly review their main conclusions. BP find that for underdoped systems the low 
frequency magnetic behavior possesses three distinct temperature regimes, with a crossover 
at T cr from non-universal mean field (MF) behavior with dynamical exponent z = 2 to z = 1 
pseudoscaling (PS) behavior. For T > T cr , u> s f ~ l/£ 2 , and the product, XQi u sf ~ u sfi 2 = 
ojq is independent of temperature, while between T cr and T*, it is u s ft; which is independent 
of temperature. NMR experiments show that above T cr , in the MF regime, uo s f and £ 2 scale 
linearly with T, i.e., u s f = A + BT, while between T* and T cr , the PS regime, it is u s f 
and which scale linearly with T, albeit with a somewhat different slope and intercept of 
u a f than that found above T cr . At T cr , the temperature dependent uniform susceptibility, 
Xo{T), possesses a maximum, while £ ~ 2, as has now been verified experimentally for the 
YBa 2 Cu30 6 . 9 2, YBa 2 Cu40 8 , and La 2 _ a; Sr :c Cu04 compounds.il Below T*, in the pseudogap 
(PG) regime, £ becomes independent of temperature while u s f, after exhibiting a minimum 
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near T*, rapidly increases as 1/T as T decreases toward T c . In addition to the change in 
the magnetic fluctuation spectrum, the PS regime is characterized by a strong temperature 
variation of the quasiparticle band structure, resulting in a FS evolution,^ This evolution 
has non-trivial consequences for the transport, as discussed in the next section. 

From a magnetic perspective, the so called optimally doped systems (e.g., YBa 2 Cu 3 6 .93 
and Lai.85Sr .i5CuO4) are a special case of the underdoped systems, in which T* is compara- 
tively close to T c . As may be seen in Table |, T cr varies rapidly with doping, and approaches 
T c as the system approaches the overdoped regime, while T*, which varies relatively weakly 
with doping possesses a shallow maximum (~ 200K in the 123 system, ~ 100K in the 214 
system) as one varies the hole content in the planes. 

From this same magnetic perspective, overdoped cuprates are defined as those for which 
T cr < T c . For these systems, then, since £ < 2, the antiferromagnetic correlations are 
comparatively weak, Xo(T) is at most weakly temperature dependent, while u s f oc £ 2 follows 
the linear in T behavior found in the underdoped systems above T cr . Examples of overdoped 
systems are the T c ~ 40K Tl 2212 system and La2_ x Sr x Cu04 for x > 0.2. 

As we shall see, the crossovers seen in the low frequency magnetic behavior possess to a 
remarkable extent their charge counterparts in magnetotransport and optical experiments. 



In this section we estimate the transport coefficients in NAFLs using a simple pertur- 
bation theory. We calculate the scattering rates for quasiparticles near the FS, assuming 
the effective interaction V e ff, Eq. (|j), and then use these to calculate the longitudinal and 
Hall conductivities. We discuss the behavior of these quantities in the relevant doping and 
temperature regimes. 

Within the Boltzmann-Bloch formalism the relaxation rate for a particle near the FS is 
generally given byS 



III. PERTURBATION THEORY APPROACH 




(6) 
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where $*; is the quasiparticle distribution function displacement obtained by solving the 
linearized Boltzmann equation. To second order in the coupling constant the scattering rate 
between points k and k' is given by 

M(k, k') = 2g 2 Im X (k - k', e' - e) [n(e' - e) + /(e')] (7) 

where e = e(k), e' = e(k'), n(e), and /(e) are the Bose and Fermi distribution functions 
respectively. Note that Eq. (|7|) is obtained by summing the Matsubara frequencies and as- 
suming a retarded scattering potential.0 No lifetime effects are included. One can, of course, 
easily add the self-energy corrections to the bare band spectrum e^, but the life-time effects 
are considerably harder to include and in principle should be obtained self-consistently. In 
the relaxation time approximation one replaces the quasiparticle distribution displace- 
ment, by ip{T){y ■ n), where tp is a function of temperature independent of k. The assump- 
tion that if does not vary greatly in the Brillouin zone is most certainly invalid for HTS,ill 
given the strong momentum dependence of the effective interaction V e ff, Eq. ([3]). However, 
it is still quite reasonable, as long as the temperature dependence of tp is the same in each of 
the hot spots, which is almost always the case. Since the inclusion of (v • n) produces only 
a geometric factor (vertex correction), which is large mostly for small momentum transfer, 
and no additional temperature dependence, we can approximate the relaxation rate by the 
scattering rate, which, in the second order Born approximation reads, 

-=/|^M(k,k'). (8) 

One can make a change of variables in the integral in Eq. (H) : 

d 2 k'^ Jde' Jy (9) 

and solve the integral over e' analytically for k near the FS, (e ~ 0), using: 

/ „,„(_,) _ f(x) = _ + , £ A_J_ = _ _ h (! + _) (10 ) 

where 
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h(x) 



X + 1 



2 ; -HI 



(id 



which has the following limiting behavior: 



h(l + x) « In 2 — |x 



1 i 

2 as+1 



X 



oo 



(12) 
(13) 



Thus the integral in Eq. (10) can be approximated by: 

1 1 



7T 



2\n n 



and the scattering rate, Eq. (||), is given by: 



7T 



(14) 



1 



ag 2 tv 



FS 



dk'M(k, k') 



(15) 



where 



M(k, k') 



T 2 



(16) 



|v| UJ kk t{uJ kk i + 7rT)' 

Note that in a typical FL one has an identical expression; however u ktk / is large, temperature 
independent and only weakly dependent on k, k'. Then for T <C uj kk i the scattering rate has 
the usual T 2 temperature dependence which ceases only when T ~ 0Jk,k'- m a realistic metal, 
however, this energy scale can be considerably reduced for certain k due to the presence 
of band singularities at these points. In a NAFL, however, u ktk i is a strong function of 
momentum transfer q, i.e., u^ y has a minimum for q = Q. Thus, the scattering rate M is 
maximized for points k and k' on the FS such that q = k — k' = Q. 

We now calculate the scattering rate for a quasiparticle near the FS as a function of T 
and Ak, where Ak is the distance to a nearby hot spot along the FS. For the commensurate 
case and a large FS the hot spots are usually distributed along the FS in such a way that 
one can write, for k — k' w Q, that (k — k' — Q) 2 w (Ak) 2 + (Ak 1 ) 2 , where Ak and Ak' are 
the displacements of k and k' from the nearby hot spots (see Fig. |l|) (the adjoint hot spots 
correspond to Ak, Ak' = 0). This basically means that there is no strong nesting at q ~ Q, 
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and the FS and its shadow FS, obtained by translating the FS by Q, intercept at a high 
angle Q ~ n/2. Then one easily finds: 

1 _ ag 2 T^T I 1 1 \ 

T k ^Vf \^UJ s f + CJ (Ak) 2 ^LU sf + llT + CU (Ak) 2 ) 

where we have assumed that the integration over Ak' in Eq. (0) can be extended to infinity. 
Eq. flT7|) can be rewritten as: 

airg 2 T 2 y/LUo 

where 



Q k = y/u af (l + (^Ak) 2 )^ sf (l + (£Ak) 2 ) + nT x 

(^(1 + (£Afc) 2 ) + ^(1 + (£Afc) 2 ) + ttt) (19) 

Note that this result is easily generalized to arbitrary Q by replacing Ak with (Ak) sin Q. 

Before proceeding with the transport coefficients it is instructive to calculate the scat- 
tering rates for different values of Ak. Eq. ( |19"1) suggests that there are primarily two regions 
of interest: we call a region of momenta k hot if 

Ak < 1/f (20) 

and we call a region cold if 

Ak > l/£ (21) 

Clearly the relaxation times in these two regions are different. In the hot region one has 



Qk ~ s/uTfy/Usf + (y/Uhf + ^uj sf + ixT) (22) 

The fits to NMR experiments of Barzykin and Pines,il show that in all HTS one can assume 
for practically all temperatures of interest that, to a good approximation, 

u ef < ttT (23) 
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and hence in the hot region we have 

Tk 4:V F 

On the other hand, in the limit £ ^> 1, in the cold region one has 



(24) 



Q k « y/u (Ak) VttT + cu (AA;) 2 (^o(AA ; ) 2 + ^vrT + cu (A£;) 2 ). (25) 
Here we distinguish two different cases, depending on the value of Ak, yielding: 

\ nT»u (Akr (26) 

7 = ^(fifc ^ « . (AA;) 2 (27) 

Obviously, the scattering rates depend on the particular temperature regime for the 
magnetic fluctuations, through the temperature dependence of u a f and uq. In the z = 2 MF 
regime ojq = w s /^ 2 = const and uj s f = A + BT, and therefore it is easy to verify (see Eq. 
(|1~9|)) that the behavior of the scattering rate for any Ak can be obtained from the Ak = 
case by replacing A in uo s f with A + ^(AA;) 2 . In other words, every point on the FS has 
a similar shape of l/r^; the hot region is enlarged with increasing temperature, as can be 
verified from Eq. (|2~T|). 

In the z = 1 PS regime the situation is somewhat different, due to the temperature 
dependence of u , u ~ 1/T. First, from Eq. (p4|) it is clear that the scattering rate due 
to spin-fluctuation scattering is roughly temperature independent in the hot region. The 



scattering rate in Eq. ( pq) , valid for £Ak ynT/cu s f, is applicable in this case as well. 
Finally, far away from hot spots, Eq. (0), the scattering rate is proportional to T 3 . Note 
that the condition uj Q (Ak) 2 ^> 7rT is satisfied for significant portions of the FS at very low 
temperatures. The scattering rates as functions of temperature for several different values 
of Ak, in MF and PS regimes, are given in Figs. [2]a and |3]a. 

In the PG case the situation is even more complicated and somewhat unclear, due to 
the FS evolution mentioned in Sec. [n[ Here Ak does not correspond to a fixed point on the 
FS, since the FS itself migrates with decreasing temperature. Nevertheless, provided that 
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the FS is still relatively large, the scattering rates are nominally the same as in the PS case, 
since u = k> a /£c ~ V^> where £ c ~ const is the correlation length at T c , and u s f ~ 1/T. 
However, the scattering rate at a particular point in the Brillouin zone is non-trivial and we 
do not present it here. 

So far we have neglected the fact that in deriving Eq. (|T7p the integral over Ak' was 
performed from (—00,00). This is quite reasonable provided that the system is far from the 
PG regime. However, as mentioned in Section [H], as the temperature is lowered toward T* the 
quasiparticle spectral weight is removed from the hot spots due to strong spin-fermion vertex 
corrections and at T ~ T* the system actually starts to lose parts of the FS. Hence, one 
must assume a finite limit of integration over Ak', yielding a somewhat more cumbersome 
form of Eq. (|T7|) which we do not quote here. However, it is straightforward to show that 
although this does not affect the cold regions very much, the scattering rate in the hot region 
is considerably modified, as shown in Figs. |2|b and |3|b. Although somewhat surprising, this 
result is natural in view of the FS modulation: the hot spots migrate with decreasing 
temperature and the effective size of the FS is reduced, thereby reducing the scattering of 
quasiparticles in the magnetic channel. More importantly, since the FS is small, i.e., Ak max 
is very small, the scattering is anomalous around most of the FS, as can easily be verified 
from Eq. (|19|). We return to this point shortly. 

The low field longitudinal and Hall conductivities are calculated according to: 

a xx = -2e 2 ]T(v k -n) 2 T fc ^ (28) 



a xy = -2e 3 J2 ( v k ' n -L r k) v fc x B • V(v k ■ nr k ) (^j^J 



(29) 



where n and are unity vectors perpendicular to the applied electrical field E, and B is 
the applied transverse magnetic field. 

We start with the resistivity. At relatively low temperatures, T <C Ef, where Ef is the 
Fermi energy, one can perform the change of variables @ and integrate over e. Assuming 
that the FS is four- fold symmetric (tetragonal symmetry), the conductivity is given by: 
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e 2 



a ™ = ^ J FS dkTkV f (30) 

Assuming that the Fermi velocity does not vary appreciably around the FS, we find that 
the conductivity is roughly proportional to the average of the scattering time fll8|) around 
the FS. From Eq. ( |18|) it is clear that the main contribution to the conductivity comes from 
the cold regions of the FS where is the largest, while the hot spots (Ak ~ 0) contribute 



relatively little. On averaging Eq. ([OJ) we obtain a rather cumbersome formula for the 



resistivity, which, to within logarithmic corrections, can be approximated by: 

Pxx v f T + T(Ak max ) 2 [6L) 
where T is a crossover temperature given by 

T ps u (Ak max ) 2 /2ii = (^Ak max ) 2 uj s f/2n. (32) 

If we define cold spots as these points on the FS which are maximally distant from a hot spot, 
then Ak max is the distance between a cold spot and the nearest hot spot. As may be seen in 
Fig. |, the size of the Ak max is related to the departure of the FS from the magnetic Brillouin 
zone boundary. For this FS and any large FS, such as is found by ARPES measurements 
for YBa2Cu307 and related compounds, Ak 

max ~ 1- On the other hand, for a profoundly 
small doping level, or in the presence of substantive incommensuration (which increases the 
number of hot spots and acts to reduce significantly the distance between hot and cold 
spots), one finds Ak max < 0.3. 

It is evident from Eq. (|3T| ) that in the absence of any evolution of the Fermi surface with 
temperature, p xx will be linear in temperature as long as T is appreciably greater than To. 
In the MF regime, T is independent of temperature, while in the PS regime it increases 
as £(~ 1/T) as T decreases. If Ak max and £(T*) are small enough, then T < T*, and 
the resistivity will be linear in T down to T*. The exact condition for this extended linear 
in T behavior reads (Ak max ) 2 ! £(T*) < (T*/100)if(c/53meV), where c is the pseudoscaling 
spin velocity, uj s f/£. Under these circumstances, any change in p xx at T cr will be marginal, 



since it can only arise from changes in the hot region contribution to p xx . The linearity in 
T continues above T cr until T ~ uj , an experimentally almost unreachable temperature. 
Moreover, we have assumed that the only relevant scattering is in the magnetic channel. 
Strictly speaking, in addition to the contribution proportional to \ defined above one must 
add the FL contribution, Xfl, present in all FLs, to the effective interaction V e ff The latter 
is likely to dominate at these temperatures. Without this term at T > uoq one then has 
p xx oc T 1 / 2 , provided lu s j oc T. On the other hand, from Eq. fll9l) it is readily shown that 
if u s f ~ const then p xx oc T 2 . We shall return to this point when discussing overdoped 
samples. 

What happens at T*? As mentioned in Sec. |J, i n the vicinity of T* the quasiparticle 
band structure begins to acquire a considerable temperature dependence, implying a FS 
evolution as well. Recently, Chubukov, Morr and Shakhnovich@ studied the FS evolution as 
a function of the coupling strength g at T = and found that with increasing spin-fermion 
interaction the FS acquires features of a preformed spin-density-wave state, as depicted in 
Fig. ^. Alternatively one can show@ that a similar evolution occurs at fixed g, but as a 
function of temperature, associated with a relatively strong temperature dependence of the 
spin-fluctuation damping.^ As a result, the FS appears to get closer to the magnetic Bril- 
louin zone boundary, i.e., the effective size of the FS decreases, with decreasing temperature. 
Hence when calculating transport coefficients one should assume that Ak max acquires a tem- 
perature dependence, the precise form of which is not well known at present. Below T* one 
finds that u retains its 1/T overall temperature dependence although it changes its slope: 
£ becomes constant, while u s f increases as T _1 . Whether one gets a crossover at T in the 
PG regime will then depend sensitively on the interplay between Ak max (T) and u> s f(T). If 
a crossover to p xx cx T 2 /T occurs, the resistivity will then decrease faster than T 2 , due to 
the temperature dependence of T . 

Let us examine the crossover in the vicinity of T* in more detail. The above mentioned 
FS evolution usually closes channels for scattering as the temperature is lowered below T*. 
New hot spots are developed and a smaller portion of the FS is now cold. At the same 
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time uq is strongly temperature dependent in both PS and PG regimes. Assuming that 
T ^> Tq(T) above T* the reduction in size of the cold region implies that a large part of 
the FS has a scattering rate such as those shown in Fig. |3| for small values of Ak, where 
1/t is roughly independent of T. At a certain temperature TJ the FS is effectively so small, 
that the quasiparticles around most of the FS exhibit precisely these (almost constant in T) 
scattering rates; the resistivity thus acquires a finite negative curvature, i.e., d 2 p xx /dT 2 < 0. 
A signature of this should be seen in the T variation of the scattering rate at a hot spot 
Ak ~ 0. In the next Section we shall see that this is indeed the case. It is important 
to realize that this can occur only when ujq is strongly temperature dependent; hence TJ 
is never larger than T cr , where uq = const, although in principle it can be comparable to 
T cr . A rapid temperature dependence of u Q eventually leads to a Fermi liquid like behavior 
at T TJ, regardless of the precise size of the FS. Note that perturbation theory is only 
marginally applicable in the PG regime, since one must include the strong vertex corrections 
at or near hot spots, although the fact that the conductivity is dominated by regions away 
from hot spots, where the above corrections are not as pronounced, allows us to make this 
qualitative argument. We note that the present theory is somewhat equivalent to the work 
of Quader et al,0 although in their case the scattering channels which open up at TJ are 
attributed to the multiband structure in bilayer cuprates and the proximity of one of the 
bands to the Fermi level. 

Overdoped materials exhibit only z = 2 scaling and in general u s f and £ show very 
little T dependence. Quite generally in overdoped samples u a f = A + BT displays a much 
larger value of A, associated with the strong fermionic damping of the spin fluctuations, 
found at high doping levels. Hence, even at a hot spot (Ak = 0) the scattering rate I/t* is 
only weakly anomalous, i.e., 1/t& ~ T or even T 2 . Therefore in the overdoped samples one 
encounters a crossover from quadratic in T behavior of the resistivity to linear in T behavior 
above To, even though the amount of the anomaly in the hot regions is weak. However, there 
are extreme cases, such as the overdoped samples of 214 (x > 0.25), where the anomalous 
behavior is so weak that the resistivity is never linear in T. Moreover, in overdoped materials 
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both the large values of u s f and the small values of £ imply that the anisotropy of scattering 
rates in overdoped materials is seriously reduced (see Eq. fll9|), producing only a weakly 
temperature dependent Hall coefficient. 

We now turn to the Hall effect. Here the situation is, in principle, a lot trickier, since the 
Hall coefficient is a non-trivial function of temperature, due to the strong anisotropy of the 
mean free path, £ = r^Vk, as readily observed from the definition of a xy (Eq. (p9|)). However, 
in practice complications occur only if I has a very peculiar periodicity as k goes around the 
FS. For the FSs seen in ARPES measurements on a number of samples near optimal doping 
one can safely assume: 

e 3 B f 

a -y = ^ J FS dkfk (33) 

where we have performed an integration over e in Eq. (^9]) in the usual manner. Again, 
assuming that Vf does not vary appreciably around the FS, one concludes that a xy is roughly 
proportional to the average of r|. Clearly, the leading contribution to this quantity comes 
from the regions of the FS where Ak ~ Ak max (cold regions); the hot regions contribute 
very little. On averaging r| over Ak (Eq. (|TT?D), we make the same approximations as before, 
i.e., uj s f <C tiT and Afc£ ^> 1 in the cold regions of interest. The result for a xy is again 
cumbersome, but can be simplified by expanding a xy T 4 in T, leading to a relatively simple 
formula: 



, 5 3 R„2 (Ah \3 



e 6 Bvj (Ak n 



(3T 2 + 6T T + T 2 /V2) (34) 



where T is defined in Eq. (p2|). Eq. (0) is valid over an extended temperature range. 
Clearly, through the variations in T and Ak max , it leads to a rather complex morphology 
of the Hall effect. 

Our result, Eq. (|34|), suggests that there are three temperature regimes for the Hall effect 
in NAFLs: In the low temperature regime a xy oc Tq/T 4 , in the high temperature regime 
a xy oc 1/T 2 , while at the intermediate temperatures a xy oc T /T 3 . The exact values of the 
crossover temperatures depend greatly on the details of the band structure. In what follows 
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we shall assume, for the purpose of simplicity, that the crossover between the low and inter- 
mediate temperature regimes occurs at T and that the crossover between intermediate and 
high temperatures occurs at 27rT , where T , given by Eq. fl32|), is the crossover temperature 
observed in the resistivity; the estimate of other quantities, such as the Hall coefficient, is 
then rather convenient. Our results for the transport coefficients are given in Table [TT]. Note 
that for many systems away from PG regime one finds cot Oh ~ T 2 , due to the relatively 
large value of the second crossover temperature T2. 

Our discussion of the resistivity in underdoped cuprates is easily extended to a xy \ in the 
MF regime To is constant and hence 

<? xy ~ T' 3 (35a) 
cot H ~ T 2 (35b) 

p xy ~ f/T (35c) 

for temperatures well below u (Ak max ) 2 . At somewhat higher temperatures the second 
crossover in a xy is approached so that both p xy and cot Oh deviate from this behavior. In the 
PS regime a xy oc T /T 3 , decreases faster than I/T~ 3 , since T decreases with temperature. 
As T approaches T cr one finds a xy ~ T -3 even in this regime. Note that in this case the Hall 
coefficient varies faster than 1/T. Finally, due to the variation of T , the second crossover 
can be approached, in which case p xy becomes only weakly temperature dependent. We 
emphasize once more that in both MF and PS regimes the results depend a great deal on 
details of the band structure. Moreover, the presence of impurities makes a non-trivial effect 
on all of the quantities of interest. 

As before the PG regime is even more complex due to the very rapidly varying u 
and the FS evolution, which alters Ak max . Although it is not clear what the temperature 
dependence of T is in this case, due to the fact that at present one does not know the 
temperature variation of the band structure in detail, we can still make certain arguments. 
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Near T* the value of To changes rapidly and one can easily switch between several regimes, 
due to the sensitivity of the Hall effect to the variation in the band structure. For example, 
one can jump from the first (low T) regime to the third, without noticing a substantial 
difference in, e.g., the Hall coefficient, since in both cases the quantity is roughly constant, 
while a xx and a xy both experience non-trivial changes. Finally, at very low temperatures 
one finds p xy approximately constant and cot G# varying faster than T 2 . 

The NAFL model suggests a very complex morphology of the transport coefficients in 
HTSs. We explore the extent to which this morphology is experimentally supported in the 
next Section. 

IV. ANALYSIS OF TRANSPORT MEASUREMENTS 

In this section we review briefly the in-plane normal state transport properties of the 
cuprate superconductors, with particular attention to the extent to which the doping and 
temperature dependence of the planar resistivity, p xx , and the Hall conductivity, a xy , reflect 
the magnetic crossovers discussed in Sec. [TT] and the perturbation theoretic estimates of 
scattering rates given in Sec. |TXT| . We then use a phenomenological model which incorporates 
the expected highly anisotropic quasiparticle mean free path to deduce from experiments 
the detailed behavior of the quasiparticle scattering rates in both the hot and cold regions 
of the Fermi surface. We focus on the results obtained for the bilayer 123 system (of 
which YBa2Cu40s, whose behavior corresponds closely to that found in YBa2Cu306.68, is 
an "honorary" member), the single layer L^-xSr^CuC^ (214) system, and as representative 
of overdoped materials, the 15K single layer Tl 2201 material. 

It is often stated that optimally doped materials, such as Lai.ssSro.isCuCU and 
YBa2Cu3 O6.93 possess a resistivity which is linear in T over a wide temperature region 
which extends down to T c , However, close examination of single crystal data shows that as 
T approaches T*, departures from linearity occur in these materials, departures which are 
in fact characteristic of all underdoped materials, but which become more pronounced as 
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one reduces the doping level below the optimum level. Quite generally, p xx in underdoped 
cuprates exhibits a drop (below linear behavior) at a temperature, TJ, which for many sys- 
tems is not far from the temperature, T*, at which the low frequency magnetic behavior 
crosses over from the PS to the PG regime. On the other hand in overdoped materials, p xx 
displays an upturn (from linear behavior) with decreasing T, an upturn which we attribute 
to the comparatively weak anomalous scattering in the vicinity of the hot spots we have 
discussed in Sec. Thus one encounters qualitatively different departures from linearity, 
depending on whether one is describing a system which is overdoped or underdoped (using 
the magnetic classification proposed by Barzykin and Pines). Some of the experimental 
results which have led us to this conclusion are given in Fig. [|. We note that the transport 
results presented in Fig. ^|a provide additional support for the assignment, based on analysis 
of NMR data,il of "optimally-doped" YBa 2 Cu 3 6 .92 material to the underdoped sector. 

To what extent does the crossover at XJ correspond to the magnetic crossover at T*. 
As may be seen in Fig. |5|a and Fig. |6], for YE^C^O^ and YBa2Cu30"6.63 and 124, within 
the present experimental uncertainties, the two crossover temperatures agree. However, the 
perturbation theory analysis of Sec. |TJ suggests that this need not be the case in general, 
since if the FS is sufficiently far from the magnetic Brillouin zone boundary, TJ can be 
considerably higher than T*; this appears to be the case for the underdoped members of 
the 214 family (see Fig. [5]b for Lai. 9 Sr iCu0 4 and Table |). On the other hand, to our 
knowledge there are no samples for which TJ is appreciably lower than T*. 



We have called attention, in Sec. [TTJ, to the possibility of a crossover to almost Fermi 
liquid like behavior for p xx . This occurs when T < T (T) (see Eq. ([32])). This crossover 
is found in most, if not all, underdoped materials, at sufficiently low temperatures, where 
p xx displays a finite positive curvature above linearity, p xx = p xx (0) + /3T b , with b larger 
than 1. Indeed, as may be seen in Fig. |8|, reasonable fits to the data may be obtained with 
b ~ 2. Note that here p xx (0) is the actual residual resistivity, due presumably to disorder in 
a sample. 

In the overdoped materials one also observes a crossover from p = p xx (0) + f3T 2 to linear 
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in T behavior at higher temperatures, as already illustrated by the Tl 2201 result in Fig. |5|a. 
We emphasize once again that the physics of this crossover is different than in underdoped 
materials, as it arises from the fact that the the anomaly of scattering at hot spots is rather 
weak in these compounds, rather than in the pseudogap effects found in the underdoped 
cuprates. It is important to notice, however, that there is no second crossover, i.e., TJ has 
not been observed. Moreover, in some overdoped materials, such as La2_ x .Sr x .Cu04 at the 
doping level x = 0.3, the resistivity is never linear in T, i.e., p ~ T 3 / 2 over the entire T 
range S 

The morphology of the resistivity is still more complex in materials with CuO chains. 
As may be seen in Fig. |7|, where we show how p aa and pbb deviate from linearity in the 124 
compound, chains lead to both a different magnitude and a different temperature dependence 
for p aa and pbb (see Ref. |37|). Qualitatively similar results are found for optimally doped 123,11 
although the level of a-b plane anisotropy is considerably less prominent. On comparing Figs. 
H and |7| it may be seen that the temperature dependence of p aa and p^ in 124 resembles 
that found in the underdoped and overdoped cuprates respectively. A detailed explanation 
of this unique behavior in terms of NAFL theory is given in Section |V|. 

We consider next the Hall effect in a transverse magnetic field. In ordinary FLs one 
finds a Hall conductivity which is roughly proportional to the square of the longitudinal 
conductivity, io xy ~ & xx cr yy ), and as a consequence the Hall resistivity, p xy = o xy j \o xx a yy — 
@xy@yx) ~ @xy xxVyy, is only weakly temperature dependent. However, this is not the case 
in HTSsJI where typically the Hall resistivity decreases sharply with increasing temperature. 
As discussed in the Introduction, it has been suggested (in the context of the RVB ground 
state) that the scattering processes which involve p xx and p xy are intrinsically different, 
and that the cotangent of the Hall angle, cot 0^ = p xx /p xy is universal and proportional 
to T 2 (Ref. [T0| ) . Much of the experimentally available data focus on the Hall coefficient, 
Rh = Pxy/B, where B is the applied magnetic field, although the above proposal has led a 
number of authors to plot cot 0#. However, we find that an examination of the behavior 
of a xy (T), a quantity more directly related to theoretical calculations, provides more insight 
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into the underlying physical origin of the measured anomalous Hall effect behavior. 

Since p xx typically varies linearly with T for T > T , a xy must decrease faster than T~ 2 
in order to bring about the temperature dependent behavior of Rh- Examination of the 
experimental data presented in Figs. [5] and [K] shows that this is indeed the case at high 
temperatures, where a xy is seen to vary as T~ 3 for both the 123 and 214 systems. However, 
experiment also shows deviations from this high temperature behavior in many systems. 
From Eq. fl34|) it is clear that if T is temperature independent, a xy should indeed vary as 
T -3 . However, in PS regime this is not the case, since To ~ 1/T and hence the result in Fig. 
H is to some extent in contradiction with our result. However, the inset of the Figure shows 
clearly that in YBa 2 Cu30 6 . 63 a xy obeys T~ 4 quite well up to temperatures T ~ 200K, where 
the first deviations from T -4 behavior occur. In fact, only well above this temperature does 
<j xy become proportional to T -3 , in agreement with our predictions. One can observe similar 
deviations at other doping levels in this system. Note also that a xy ~ T -3 is obtained also 
well above 27rT , for T ~ 1/T (see Table [TJ). We draw two important conclusions. First, 
<j xy has a crossover from T -4 to T -3 behavior, the details of which depend not only on 
whether one is above T , but also on the temperature regimes considered. Second, and more 
importantly, the "disagreement" in the experimental data in Fig. [5] and its inset show clearly 
that one should not use simple minded power law fits, as is often done in the literature, but 
more complex formulae like Eq. ([34]), which include all necessary crossovers. In addition, 
we note that one must be very careful when fitting to the Hall effect data, since impurities 
play an important role, and their presence can lead to large positive intercepts of l/u xy as 
a function of T 3 , such as that seen in Fig. [TO. Moreover, as mentioned above, the chain 



contribution in various materials is non-trivial and may obscure a crossover in the in-plane 
Hall effect. Finally, we remind the reader that Eq. (0) is obtained in the limit u s f <C T. 
This limit is not well satisfied in highly overdoped materials where a somewhat different 
sequence of crossovers may occur. 

Further insight into the role played by anisotropic quasiparticle scattering may be ob- 
tained using a phenomenological model for the variation in the mean free path as one goes 
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around the FS. For simplicity, we consider a representative FS and anisotropic quasiparticle 
behavior which exhibits a four-fold symmetry. We assume a cylindrical FS with approxi- 
mately constant kf, but with an anisotropic effective mass m^. We parametrize the mean 
free path (MFP) by assuming that there are hot regions on the FS at 90° with respect to 
each other, and choose as a representative FS that depicted in Fig. |I|. In the extended zone 
scheme this FS is approximately a large distorted circle, centered around (vr,7r), with hot 
spots near k = ±(7r, 0) and k = ±(0, 7r). Then the MFP around the FS is given by: 



1 1 

£hot 



1 + a cos(46>) 



(36) 



1 + a 

where the anisotropy parameter, a, varies between and 1 (a ~ 1 for a highly anisotropic 
scattering rate and a « for weak anisotropy), and £h ot is the MFP at a hot spot on the 
FS. Note that, in general, both a and £h ot are functions of temperature. For simplicity 
we shall assume a = (1 — r)/(l + r), where r = £hot/£coid = £(0) / £(ir / 4) . This particular 
choice may not be the most suitable for an effective interaction which is sharply peaked at 
Q, since for long correlation lengths, as one moves away from a hot spot l/£ decays faster 
with 6 than is assumed here. However, if the effective interaction ([I]) has incommensurate 
peaks in momentum space, then the anomalously scattered (hot) regions are large, l/£ is a 
slower function of 6, and Eq. ( |36"|) is quite reasonable approximation. Therefore, one should 
consider £h ot and £ co id here in a broader sense, recognizing that under certain circumstances 
it may not be possible to identify them with the results quoted in the previous section. 

Starting from the expressions, Eqs. (|28]) and (|29|) , for o xx and a xy , at T <C t one can 
perform a change of variables and integrate over e. The result for <j xx is: 

e 2 k t 



Or 



2tt 2 



d6£ cos 2 6 (37) 



where £ = r(6)kf /m(6) and 8 is the angle between the electrical field and k(#). On using 
the parameterization (pB]), we find: 

ne 2 r hot ne 2 r co i d ^/r 

Oxx — 1= — [<5o) 

m hot y/r m co i d 
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where n = k'j/2'K. In similar fashion the Hall conductivity is found to be: 

a xy = (^^j J d6l cos 9 (J^i sin ej (39) 

which leads to: 

/ e 3 B\ 2 1 + r 1 + r 

^ = J = °^ B Zhot—. (40) 

Eqs. ( |3lf ) and ([K]) may be combined to yield simple expressions for the Hall coefficient and 
Hall angle: 

a xy 1 1+r , s 

RH "^fB = Tn-^fi (41) 

and 

ootfe = f« = _^*2L JlL (42) 
CTxj, eBr hot 1 + r 

Note that the above expressions depend only on the two temperature dependent parameters 
r and ihot- These two quantities can be directly probed by measuring the Hall coefficient 
( p| ) and the resistivity p xx ~ o~~ x l , since Rh depends only on r and the product p xx Rh yields 



<-hot- 

In the case of high anisotropy, i.e., r < 1, the above quantities take an especially simple 
form: 



2 / ^coldThot ;i \ 

° xx = ne J (43) 

V m hot m co i d 

^cTcold (AA\ 

0~xy ®~xx ^ \ ) 



1 



T cold j™ cold 



Rh = — < (45) 



" 2ne\ T hot /m hot 



2 

cot9 H = (46) 

UJ c T cold 

where ui c is the cold cyclotron frequency, u c = eB/m co id. Our result, Eq. ([P]) for the Hall 
angle, provides a natural explanation of why the superconducting cuprates with a variety of 
different behaviors for p xx (T), often display a quite similar behavior of cotO#. Note that 
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the Eq. ( |43|) is exact for any r. It is clearly consistent with our conclusion in Sec. |TTT| that 
the linear in T resistivity in NAFLs is due to a fine balance between cold and hot regions. 
In addition, with the help of Eq. fl2"T|) one easily verifies that Eq. ([4*4]) is consistent with Eqs. 



(0) and 

One can use Eqs. (P8|)-(p2D to deduce the temperature dependence of (m/r) hot and 
(m/T) co id from the experimentally measured p xx and p xy . Only one undetermined parameter 
enters into this procedure, namely n, which in some cases can be obtained experimentally 
from the high temperature limit of Rh, since in this case r — > 1, and the Hall coefficient 
is temperature independent. The quantity. Roo = Rh(T — > oo) is well defined for many 
overdoped materials,0 but has been determined for only a few underdoped materials. When- 
ever experimental data are not available we make the best estimate for R^: for example, in 
the 214 family we extrapolate Rh{T) to T = T cr and for the 124 material we assume R^, 
comparable to the result found for 

YBa 2 Cu 3 6 . 6 3 (Ref. |3~|), which has comparable magnetic properties. However it is im- 
portant to stress that as long as one has the correct order of magnitude for R^, qualitatively 
the results are virtually the same. 



We proceed to extract quasiparticle scattering rates by fitting Eqs. (|5E|)-(|4"T|) to the trans- 
port measurements considered above. Fig. [ll| shows the scattering rates for YBa2Cu06.63 in 
the cold and hot regions of the FS (in the same dimensionless units) as a function of tem- 
perature. Clearly in the hot region the scattering rate is much larger; it is approximately 
linear in T for T < TJ ~ 200K, and becomes weakly T dependent at higher temperatures. 
As noted above, for this system 7J m T*. At the same time (^/r) co ^ barely reflects the 
crossover at TJ as may be seen in Fig. [12|, although clearly both scattering rates are affected 
as the system approaches the PG regime. Note that these scattering rates closely resemble 
those shown in Fig. [|. In particular, the crossover to a constant scattering rate of (m/r)^ 
is strikingly similar to the calculated one for Ak ~ 0, shown in Fig. f3[ 

The results for the 124 compound are given in Fig. [T3|, where we have taken p aa , which 
does not include chain contribution and exhibits an in-plane resistivity resembling that in 
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underdoped materials, and used the measured values of Rr- The results are similar to those 
found for YBa 2 Cu30 6 . 68 , although in the measured temperature regime (T < 300K) the 
scattering rate at the hot spot does not quite saturate. However, we find that, as T increases, 
(m i V) cold becomes approximately proportional to T 2 in this material; on extrapolating this 
behavior of (m/r)^^ and the experimentally obtained resistivity, we find the same behavior 
of (m/T)hot, albeit with a somewhat higher value of XJ ~ 250K. 



Finally in Fig. |TJ we show the scattering rates in the cold and hot regions for YBa2Cus07: 
here the scattering rate in the hot region is once more seen to saturate, but at a much lower 
temperature corresponding to TJ « 120 — 130K. Both this value and that found in 124 
and YBa2Cu3 06.68 are very close to T*, found by Barzykin and Pines for these materials 
(see Table |). Although it is not easy to verify directly in NMR experiments, because both 
T cr and T* are close to T c for YBa2Cu307, this result provides a further indication of the 
presence of the PG phase even in optimally doped 123 materials. 

A similar analysis can be performed in single CuO layer materials: two panels of Fig. 
IB] show the scattering rates at doping levels x = 0.10, x = 0.15 in 214 compound. Much 



like YBa 2 Cu 3 6 . 68, one observes the saturation of the scattering rate in the hot region 
and it is only for a particular anisotropy level that the resistivity in the optimally doped 
214 material is approximately linear in T below TJ. In fact, as pointed out in Ref. a 
crossover in resistivity exists even in this material, but is so minor that one usually assumes 
that p xx oc T at all temperatures. As mentioned before, in the 214 family 7J is considerably 
larger than T*. 

Finally we consider the overdoped single layer materials: Fig. [16| shows the scattering 
rates for the overdoped 15K Tl 2201 compound. First, one immediately notices that the 
anisotropy of scattering is not nearly as pronounced as in previous graphs. Note that in this 
case Roo is very well known and thus one can consider this reduced anisotropy in scattering 
rates quantitatively as well. We remark that both scattering rates have nearly as strong T 
dependence as in the previous cases, over the entire T range. Moreover, the scattering rates 
shown here again closely resemble the calculated result for larger Ak in Fig. |2|. 
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V. NUMERICAL RESULTS 



In the previous sections we have established the relationship between the crossovers in 
the magnetic spectrum of cuprates and their transport properties analytically. However, it 
is unclear whether such an analysis, based on back-of-the-envelope calculations, can account 
for the experimental results quantitatively. In this section we present our numerical results 
for the transport coefficients in NAFLs: we show that for spin-fluctuation and band parame- 
ters chosen from NMR, INS and ARPES measurements, the NAFL model indeed reproduces 
experimentally obtained results. We compare our results to experiments performed on sev- 
eral compounds of interest, examine the sensitivity of our results to the input parameters of 
the theory, and address several important questions regarding realistic materials, including 
the strong a-b plane anisotropy observed in a number of bilayer compounds. 

The method used to obtain the results presented here has already been discussed 
elsewhere^ and hence we only review it briefly. We study the transport in NAFLs by solving 
the Boltzmann equation (BE) numerically. We obtain the displacement of the quasiparticle 
(fermionic) distribution function, 



as a function of momentum k, temperature and frequency on a fine mesh of points in the 
Brillouin zone, using a standard relaxation method. The transport coefficients are then 
obtained from 



where v = Ve^. We assume that the interaction present in the collision integral in the 
BE, is given by Eq. ([!]). The numerical method can be described as follows: we start from 
a reasonable choice of which we then use to calculate the collision integral, which in 
turn provides for a new value of The iteration procedure is stopped when the difference 
between two subsequent values of $fc is smaller than a given tolerance. 




(47) 




(48) 
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In our earlier numerical workjHMll we demonstrated the viability of a NAFL description 
of transport in the cuprates. By using spin-fluctuation spectra seen in NMR experiments, we 
were able to obtain both qualitative and quantitative agreement with the analysis presented 
in Sec. [TTT] and Sec. [TV]. However, for the band parameters, doping level, and spin-fluctuation 
spectrum we assumed for optimally-doped YBa 2 Cu 3 7 , the crossover temperature, T a , was 
high compared to that seen experimentally, while p xx exhibited FL behavior up to temper- 
atures ~ 200K. Because T Q was large, on extrapolating our linear in T results for p xx to 
zero temperature, we found a large negative intercept. Moreover, although the calculated 
Hall coefficient decreased with increasing temperature, it displayed a far weaker tempera- 
ture dependence than that seen experimentally. Our results thus resembled those found for 
overdoped cuprates, rather than an optimally doped sample. 

On the basis of the analysis presented in |T|, we can identify several potential causes 
for these discrepancies. First, the actual value of u = ujsfC, 2 in optimally doped 123 
is < 60meV,ii rather than the larger value, c^o = 76meV, adopted in these calculations. 
Second, the spin fluctuation spectrum was assumed to possess a commensurate peak, rather 
than the four incommensurate peaks which a recent analysis of NMR and INS results suggests 
might provide a better fit to experiments on the system.il Third, the doping level was taken 
to be n ~ 0.25%, while the LDA calculations of bonding and antibonding in bilayer materials 
discussed below suggest a considerably lower doping level might be more appropriate. Each 
of these effects acts to bring about either a lower value of T , an increase in the linear in T 
regime of p xx , and/or an increased temperature dependence of Rh- To explore their combined 
influence, we carry out a numerical solution of the BE for a representative "standard" 
123 material for which n = 20%, ojq = 60meV, and u>sf(K) = 60 + 0.6T(K); the band 
spectrum was assumed to be (t = 0.2eV; t' = — 0.35t), which is slightly different from the 
unrenormalized LDA spectrum,@ while the incommensurate peaks in the spin spectrum are 
assumed to lie along the diagonals, at Qj = Q + (±0.l7r, ±0.l7r), as proposed by Zha et al.il 

In Fig. [T7] we compare our calculated resistivity as a function of temperature with the 
experimental results obtained for optimally- doped YBa2Cus07. The coupling constant g, 
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which sets the scale for the magnitude of the resistivity, has been adjusted to g = 0.48eV, in 
order to obtain the best fit. Clearly, the calculated result agrees very well with the data down 
to temperatures only slightly above T c . Note that the extrapolated value of p xx at T = 
is now quite small, as is seen in the experiment. On the other hand, the Hall coefficient, 
shown in the inset of the figure, is strongly temperature dependent, with a relative decrease 
which is approximately the same as that seen experimentally in the temperature interval, 
150-400K, although the exact shape of Rh is not the same. However, as mentioned above, 
a better knowledge of the band structure, including the chain contribution, and of the 
effective interaction V e ff for q away from Qj might well remove this problem. Moreover, a 
quick look at Eq. fl45|), as well as Eqs. ( |31~D and (|34"D indicate that the FS reduction (i.e., 
decreasing Ak max ) with the approaching PG phase at low temperatures, should provide for 
an additional temperature dependence of Rh- Our calculated Rh saturates at higher T, to 
a somewhat higher value than that seen experimentally, implying that a different value of 
the doping level should be used; nevertheless, our choice of parameters would seem quite 
reasonable (see Eq. (|4"T|). Another reason for a discrepancy in the shape of Rh lies in the 
fact that we have calculated this quantity assuming a tetragonal lattice, rather than the 
orthorhombic one found in YB^CuaOy. Recent measurements^ on twinned and detwinned 
crystals of this material show a considerable difference in Rh, indicating the importance of 
including the a-b plane anisotropy in the calculation. We return to this point at the end of 
this section. 

The doping level used in this calculation is somewhat low. If one assumes that every 
Oxygen atom doped into the 123 system extracts a single hole from CuO planes, the doping 
level of the optimally doped 12306.93 should be somewhat higher, n « 23%. On the other 
hand, simple arguments, based on LDA calculations, suggest that in bilayer materials, such 
as 123 compounds, there are two bands, bonding and antibonding, separated by twice the 
hopping matrix element between the CuO layers, t± (see Refs. |4l3||42| ). If one assumes that t± 
is weakly momentum dependent, one finds the top (antibonding) band is more heavily doped, 
while the bottom (bonding) band should be close to half filling, regardless of the Oxygen 
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content .El Since the interaction between quasiparticles is much weaker in the metallic state 
of cuprates than it is in the insulating state, the nearly half filled bonding band has a larger 
value of vj and should provide the dominant contribution to the conductivity in the bilayers 
(see Eq. (|3~ID). Hence the doping level chosen here, n = 20%, may be quite reasonable for 
123 and 124 compounds. 

To study the effect of doping on the resistivity at fixed spin fluctuation spectra we 
have calculated p X x{T) for the same values of the spin and band parameters as in our 
"standard model," but with the chemical potential adjusted to yield an n = 30% doping 
level, and compared it to the above case in Fig. |l^. Clearly, 

Pxx CGclSGS to be linear in T 

at lower temperatures, while continuing to exhibit typical FL behavior as T — > 0. In the 
same figure we also show the influence of incommensuration, by plotting p xx (T) for the 
same representative set of parameters as above, but with a commensurate spin-fluctuation 
spectrum. Although the commensurate result has a slope at higher temperatures which is 
remarkably similar to the incommensurate (representative) one (due to the fact that p xx 
is independent of u above T (see Eq. (0)), we again observe FL-like behavior at lower 
temperatures. In both cases the explanation is simple: enlarging Ak max , leads to an increase 
in T ; this can be achieved either by increasing the doping level, which increases the size 
of the FS, or by removing the extra hot spots brought about by incommensuration. Fig. 



18| demonstrates quite clearly that the behavior of the resistivity is sensitive to the choice 
of band parameters and leads us to conclude that self-consistent calculations which take 
into account changes in the quasiparticle interaction should be used in comparing with 
experiment .0 On the other hand, assuming that the magnetic (band) parameters are well 
known, our model can be used to put constraints on the allowed band (magnetic) parameter 
values. 

Next we examine the dependence of p xx and a xy on u s f = A + BT. In Fig. |19| we plot 



the resistivity as a function of temperature for four different values of A and B: we show 
the cases where A = 20K and B = 0.2, A = 0.6 and B is the same, the case where B is 0.6 
and A is the same, and the case where w s f = 60 + 0.6T. The band parameters are set to 
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t = 250meV, t' = — 0.45t and the doping level is n = 15%; the coupling constant is set to 
g = leV. In agreement with our analysis in Sec. |TTT| , a large value of A, corresponding to an 
enhanced damping of spin-fluctuations due to spin-fermion interactions, shifts the resistivity 
down; it leads to a higher value of To in Eq. fl3~T|) and only marginally anomalous scattering 
rates even in hot regions of the FS, in agreement with the result presented in Fig. 

On the contrary, increasing the value of B only changes the slope of the resistivity curves 
shown in Fig. [19|, but does not much affect the behavior with respect to the origin. This is 
to be expected, since according to Eq. fll9|) the term BT in u s f acts in conjunction with the 



7rT term in yvrT + cj s /(1 + (Ak) 2 ^ 2 ) to bring about linear in T resistivity. However, since in 
all experimentally studied samples B <^ n, the effective value of T ~ u> (Ak max ) 2 /2(n + B) 
differs only slightly from that quoted in Eq. (|32|). As noted in Sec. |TJ, for practically all 
cuprates BT is never large enough to produce the T 1//2 behavior of the resistivity discussed 
there, which is found at T ~ LU (Ak max ) 2 , provided the FS is large enough. Hence, the 
crossover at TJ, where a curvature resembling T 1//2 is found in many underdoped samples 
at T > T*, has to be attributed to a small value of Ak max resulting from the FS evolution. 
Finally we note that the Hall conductivities for all of the above cases obey a T~ 3 law quite 
well, as shown in the Inset of the figure, although due to the large value of To the Hall 
resistivity may be only weakly temperature dependent. 

We turn now to the z = 1 PS regime. As indicated in Sec. |HI| , the resistivity is linear 
in T for T > T regardless of whether or not T is temperature dependent. However, 
since £ ~ 1/T in this regime, the condition £A/c < 1 is satisfied over most of the FS at 
relatively low temperatures. In this case the resistivity due to spin-fluctuation scattering 
tends to a constant, since T ^> 2ttT . To further explore the PS regime, we solve the 
B.E. numerically for a representative member of the 2-1-4 system, [c = u s f£ — 50meV, 
£-! = 0.1 + T/1000K, g = leV, n = 15%, t = 250meV] and several values of t' (-tf ft = 0.2, 



0.3 and 0.4 respectively). Our results are shown in Fig. [20|a. In all of the cases displayed 
there we have assumed that the effective interaction has incommensurate peaks at Qj = 
Q + (±AQ, 0) and Q + (0, ±AQ), with AQ = 0.257T, in order to be consistent with neutron 
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scattering experiments performed in 214 materials at sufficiently high doping levels.tB For 
t' = — OAt the FS is relatively large, as is the corresponding T (T). As a result one observes 
a crossover to linear in T resistivity at T ~ 200K from the low T Fermi liquid like behavior. 
For t' = —0.2t the FS is always close to the Brillouin zone boundary and hence the effective 
Ak max is small, Ak max ~ 0.1. Hence the condition T = 2irT (T) is satisfied for T ~ 1000K, 
and the resistivity deviates from linearity already at T = 150. Obviously, in the intermediate 
case Pxx o^- T over an extended temperature range. We conclude that the resistivity is 
very sensitive to changes in band parameters, especially when To is a strong function of 
temperature. 



The parameters used in Fig. 20 a are close to those found for Lai. 8 5Sr ai5 Cu04 material 



in the PS regime.c§EI In Fig. [20| b we compare p xx (T), obtained for t' = — 0.33t (solid line), 
to experiment. The coupling constant has been adjusted to yield the same slope as the ex- 
perimental curve. Obviously even in this case the calculated resistivity exhibits the features 
shown in Fig. pO|a, but the fit is quite good, especially keeping in mind the fact that the 
experimental result is subject to disorder. On the other hand, when one compares cot Oh in 
these two cases (solid line in Fig. f2T|), the fit, although remaining reasonable, is not as good. 
The curvature of the calculated result (and lack thereof in the experimental result) implies 
that the cold regions of the FS are too small in the calculation, whence cot Oh ~ T 1S , 
rather than the T 2 , behavior frequently seen, although some samples of this material have 
the smaller exponent in cot Oh found here. However, even for a small value of t', the ad- 
dition of a Fermi liquid term, Xfl = Xo(T)/(l — iiruj /T(T)), to the effective interaction, 
makes the resistivity appear linear in T up to higher temperatures and the agreement with 
experiment for the Hall effect much better, as shown by the dashed line in Figs. |20] and [21]. 
In both cases we have used T(T) = (0.5 - T/2500)eV and Xo(T) = 1 + 0.4T/100states/eV, 
the values obtained from the fits to NMR measurements.^ This is as expected, since adding 
a Fermi liquid term affects predominantly the cold regions of the FS and its effect is similar 
to that of the strong band renormalization present in the PS regime. 

We next consider the issue of the a-b plane anisotropy. Although found in practically 
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all cuprates, this anisotropy is especially pronounced in materials with CuO chains, leading 
to quite different transport results along different crystallographic directions. In principle, 
these are two reasons for an anisotropic resistivity: the chains can form a conducting band, 
or the chain bands hybridize with Oxygen bands in CuO plane, producing an anisotropic 
quasiparticle effective mass. Since the chains in, e.g., YBa 2 Cu 3 7 are very close to the CuO 
planes, and the Oxygen p-orbitals are rather large, we believe that the latter effect must 
dominate. In order to simulate its effect on the in-plane resistivity we consider a 2D band 
with anisotropic hopping integrals: we alter the first term in Eq. (g) to read 

— 2t(cos k x + a cos k y ) (49) 

In Eq. ( f49"|) the anisotropy parameter a, lies between and 1; obviously, a = 1 in tetragonal 
systems. A FS for an anisotropic system with a = 0.55 and the doping level n = 10% 



is shown in Fig. |22|. Although this shape of the FS is somewhat incompatible with that 
observed in, e.g., the 124 material, one must keep in mind that strictly speaking a single 
band model, such as that considered here, is an oversimplification; we only use the present 
model in order to qualitatively explain the physics in the above material. 

In Fig. we show the resistivities for the current running along x and y directions for 



the band depicted by the FS in Fig. |22| We assume that uj s f{K) = 20 + 0.2T, which is close 
to u s f found in 124 materials. However, we retain the z = 2 scaling regime in order to make 
the comparison with previous results easier. Clearly, not only are the magnitudes of the two 
resistivities different, but their temperature dependence is also different; indeed our results 
for p xx and p yy resemble those found in the 124 material. The lower curve (p xx (T)) appears 
to have a considerable negative intercept, indicating a large value of To, as in overdoped 
materials. The upper curve, p yy (T), exhibits a crossover at T ~ 300K, similar to that 
observed experimentally at TJ in underdoped cuprates. This quite different behavior arises 
from the modified positions of the hot spots: as may be seen on examining the FS shown 
in fig. the hot spots are now located at points (l,±7r) and (n — 1,0) and symmetry 
related points on the FS, i.e., they are asymmetric with respect to the inversion of x and y 
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axis and cannot be obtained by drawing the usual magnetic Brillouin zone boundary. For 
current running along the x direction, a large segment of the FS, roughly along the (1,2/) 
direction, can be displaced in transport (the vertex function v x is large at these k vectors) 
and more importantly, these segments of the FS are far from the hot spot at (1, ±7r), so 
that it is necessary to balance the hot and cold regions, to arrive at a p xx which is linear in 
T. For current running along the y direction, the vertex v y is small almost everywhere on 
the FS (see Eq. ( f4§l) and the hot spots sit in the middle of the region of the FS where v y 
is the largest, further reducing the conductivity in this case. Effectively the FS appears to 
be small, producing the anomalous (non-linear with negative curvature) p yy {T) seen in the 
figure. Note that the Hall conductivity is proportional to T -3 , as seen in Fig. |24j. We also 
note that our result for a xy is the same to within a few percent for both current directions, 
which serves as a nice check of the quality of numerics in this case. 

We make another important comment regarding this result: unlike the case of 124 mate- 
rial, where TJ is very close to T*, so that the anomalous behavior of the resistivity occurs at 
least in the pseudoscaling (z = 1) regime if not in the pseudogap one, here we have similar 
behavior with z = 2. We again conclude that the crossover at TJ is closely related to the 
effective size of the FS and is less closely related to a particular magnetic regime.il However, 
we must keep in mind that the FS here is reduced artificially (by introducing a very large 
anisotropy). In reality, magnetic fluctuations alter the shape of the FS as the system enters 
into the pseudogap regime. Moreover, one cannot obtain a similar anomalous behavior of 
p yy using a large Fermi surface and the same spin-fluctuation parameters as those utilized 
here; and hence this result serves as a further indication that the pseudogap regime involves 
an evolution of the FS. 

VI. CONCLUSIONS 

We have performed an analysis of the in-plane longitudinal and Hall conductivities (re- 
sistivities), in terms of the NAFL model, using both perturbation theory and numerical 
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methods. From our perturbation theory results we obtained an analytical expression for 
the scattering rate as a function of temperature for an arbitrary point on the FS, as well as 
empirical expressions for p xx and p xy . The results are then used to study analytically the 
crossovers in transport in various temperature regimes. Assuming a particular geometry of 
the FS, we obtained phenomenological expressions for both p xx and p xy which we used to 
analyze the experimental results and gain further insight into the scattering rates. Finally, 
we performed a set of numerical calculations, verifying for a realistic set of input parameters 
that NAFL model yields results in agreement with experiment. 

Our analytical results display a relatively strong anisotropy of scattering rates around 
the FS, leading to a complex morphology of both p xx and p xy , in contrast to the common 
belief that in all cuprates p xx oc T and p xy oc 1/T. A close examination of the available 
experimental data shows that the universality of the temperature behavior of both p xx and 
p xy is quite marginal over larger temperature regions, i.e., only extremely clean, optimally 
doped samples of 123 show p ~ p + aT with p ~ and cot Q H ~ T 2 . Even optimally doped 
214 compounds show considerable discrepancy from this behavior in that the intercept of 
cot ® H at T = is non-negligible. All other doping levels show more or less pronounced 
deviations from the above optimal behavior at low temperatures, in agreement with the 
NAFL model. We therefore conclude that the optimally doped compounds are rather unique 
and the temperature variation of their transport coefficients is somewhat accidental. The 
calculated results for the scattering rates are found to be in detailed qualitative agreement 
with those coming from the fits to the experimental data. It is remarkable how systematically 
the scattering rates vary with temperature and doping and how precisely they reflect the 
magnetic properties in cuprates. 

Finally our numerical results show reasonable agreement with experiments for a realistic 
set of spin-fluctuation and band parameters in both z = 1 and z = 2 scaling regimes. We 
also verified that improved quantitative agreement can be achieved through the inclusion 
of strong coupling effects, which bring about a temperature dependent quasiparticle band 
structure. In addition, we showed that the NAFL model provides for a natural explanation 
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for the unusual temperature dependence of the resistivities along different directions in 124 
material. Both the numerical results and the available experimental data provide yet another 
indication that NAFL is a strong candidate for the explanation of the anomalous transport 
in cuprates. 

It is straightforward to extend the NAFL model to optical frequencies and to take into 
account the influence of impurities and we will report on these results in future publications. 
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TABLES 



TABLE I. Crossovers temperatures and spin fluctuation parameters for representative cuprates 





% 


T* 


T cr (K) 


Lai.gSro.iCu04 


33 


50 


>500 


Lai. 85 Sr .i5CuO4 


39 


85 


325 


Lai. 8 Sr . 2 CuO 4 


30 


75 


120 


Lai. 76 Sro.24Cu0 4 


25 


35 


100 


YBa 2 Cu 3 6 .63 


60 


180 


>600 


YBa 2 Cu 4 08 


80 


210 


475 


YBa 2 Cu 3 07 


90 


110 


150 



TABLE II. The transport coefficients in various temperature regimes. 







Pxx 


Pxy 


cot Qh 


T < T 


T /T 4 


T 2 /T 


const 


T 2 /T 


T < T < 2ttT 


T /T 3 


T 


T /T 


T 2 /T 


T > 2vrT 


1/T 2 


T 


const 


T 
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FIGURES 

FIG. 1. A model of a Fermi surface in cuprates (solid line) and the magnetic Brillouin zone 
boundary (dashed line). The intercepts of two lines marks the regions of the FS near (n, 0), which 
can be connected by the wavevector Qj. These regions are most strongly scattered into each other 
(see Eq. (|)). 

FIG. 2. Scattering rates, in arbitrary units, as a function of T in the mean field z = 2 regime, 
for several values of Ak, the displacement of the wavevector k from a hot spot. In both plots we 
have assumed uj s j = 0.2T, ujq = 600K. The curves correspond to (top to bottom) Ak = 0, 0.1, 
0.25, 0.5, 1 and 2. Panel (a) shows the case of a quite large FS, Ak max = 2, and panel (b) that of 
a typical small one, Ak max = 0.5 

FIG. 3. Scattering rates, in arbitrary units, as a function of T in the pseudoscaling z = 1 
regime, for several values of Ak, the displacement of the wavevector k from a hot spot. In both 
plots we have assumed uj s f = 0.2T, ujq = 1000K/T and the curves correspond to (top to bottom) 
Ak = 0, 0.1, 0.25, 0.5, 1 and 2. Panel (a) shows the case of a large FS, Ak max = 2, and panel (b) 
that of a small one Ak max = 0.5 

FIG. 4. The FS as a function of coupling constant g in a NAFL. The figure shows only one 
quarter of the Brillouin zone. The curves are obtained using the same formalism and the same 
values of spin fluctuation and band parameters as in Ref. 35. The dashed-dotted, solid and dashed 
lines correspond to increasing value of g. As shown in Ref. 29, these curves also correspond to high, 
intermediate and low temperatures respectively. Note that the effect of the interaction is to move 
the FS towards the magnetic Brillouin zone boundary (here depicted by the dotted line). Hence 
the interaction effectively reduces Ak max . 
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FIG. 5. The measured resistivity as a function of temperature, showing the deviation from 
linearity in T in underdoped, overdoped and optimally doped cuprates. Panel (a) shows the 
results from Refs. 44, 45 and 46, obtained in (top to bottom) 15K Tl 2201, YBa2Cu307, and 
YBa2Cu30"6.63 compounds. The dashed line is a guide to the eye and the arrow marks the crossover 
from pseudoscaling to pseudogap behavior at T* in the underdoped and optimally doped material. 
Panel (b) shows La2- a; Sr, r Cu04 at three doping levels (top to bottom) x = 0.22, 0.15 and 0.10 
(Ref. 36). The quantity plotted is (p xx (T) — po)/aT, where po and a are obtained by fitting the 
high T, linear part of the resistivity. 

FIG. 6. The resistivity p aa , for current running along the a crystallographic axis, as a function 
of T in YBa2Cu40g material (Ref. 37). The quantity plotted is {p a a{T) — pps)/ctpsT, where pps 
and aps are obtained by fitting pp$ + otpsT in the PS temperature regime (200K< T < 500K; 
see Table I). The arrows mark the crossovers at T* and T cr . Note that p aa retains its linearity in 
T above T cr with only minor change in slope. 

FIG. 7. The reduced resistivity, (p aa (T) — p )/aT, in YP^C^Og (Ref. 37) a and po are 
determined as in Fig. [|. The bottom (top) set of data shows p aa (pbb), where p aa (pbb) is obtained 
with current running along a (b) crystallographic directions. Inset: p aa and pbb in YBa2Cu408 
material. Note that at higher T p aa ~ 3pbb- 

FIG. 8. Pxx{T) in underdoped cuprates at low temperatures. The three sets of data are 
shifted by 1 for clarity and correspond to (top to bottom) YBa2Cu306.68> YBa2Cu408 and 
Lai.gSro.iCu04 (Refs. 36, 37 and 46). The dashed lines are guides to the eye. The quantity 
plotted is PxxiT) Pxx (0))//^-^ > where PxxiS^) <md /? are obtained by fitting p X x — Pxx 

(0) + (3T 2 

at low T. 
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FIG. 9. The inverse Hall conductivity as a function of temperature in 123 and YBa2Cu40s 
compounds. The data correspond to the measured results obtained in (top to bottom) YBa2Cu40s, 
YBa2Cu30"6.63 an d YBa2Cu30y (Refs. 37, 45 and 46). Note that at sufficiently high temperature 
all curves satisfy T 3 quite well. Inset: the result for YBa2Cu306.63> showing that a xy oc T -4 at 
low T. 

FIG. 10. The inverse Hall conductivity for the 214 compound. The data are taken from Ref. 
36 and correspond to doping levels (top to the bottom) x = 0.10, 0.15, and 0.22. The curves show 
large positive intercepts, suggesting a considerable amount of disorder and/or a smaller FS in this 
system. 

FIG. 11. The hot (top) and cold (bottom) scattering rates in YBa2Cu30g.68 i n dimensionless 
units, obtained by fitting the results of Ito et al 46 for p xx and p xy , shown in the previous figures, to 
the phenomenological formulae (^) and (^) (see text). The rates closely resemble the calculated 
scattering rates shown in Fig. ||[ The arrow marks the crossover at T*. 

FIG. 12. The scattering rate in the cold regions of YBa2Cu306.68 obtained from experimental 
results of Ito et al. 46 Note that, to within an additive constant, this scattering rate is the same 
as that found in our perturbation theoretic calculations. The arrow marks the crossover at T*, 
discussed in the text. 

FIG. 13. The hot (top) and cold (bottom) scattering rates for YBa2Cu40s obtained from 
the experimental results of Bucher et al. 37 There scattering rates resemble closely those found in 
YBa2Cu306.68) as might be expected from the similarities in their magnetic behavior. The arrow 
marks the crossover cit T% . 
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FIG. 14. The hot (top) and cold (bottom) scattering rates for YBa2Cu3C>7 obtained from the 
experimental results of Ginsberg et al. 45 Note that l/r^ot shows a crossover from the Ak = result 
in Fig. H to the same result in Fig. [2|, which should occur near T cr , marked in the figure by an 
arrow, l/r^ot weakly increases with temperature for T <C T cr (see Fig. 0). T* is the crossover 



temperature equivalent to those shown in Figs. 11 and 13. Note that 1/t co ^ is approximately 



quadratic in temperature, in agreement with our perturbation theory calculations. 

FIG. 15. The hot (top) and cold (bottom) scattering rates for two 214 compounds obtained 
from the experimental results of Hwang et al. 36 The top (bottom) panel shows the results for 
x = 0.10 (x = 0.15) doping level. In the bottom panel the crossover temperature T* is clearly 
observed in l/r^ ot , as expected in the NAFL model. 

FIG. 16. The hot (top) and cold (bottom) scattering rates for a 15K 2201 Tl compound ob- 
tained from the experimental results of Mackenzie et al. 44 Note that the anisotropy of scattering is 
considerably reduced in over doped materials, due to the smaller correlation length and the strong 
fermionic damping, manifested in a larger value of uj s f (see text). 

FIG. 17. A comparison of our calculated result for p xx {T) with the experimental result of 
Carrington et al. 20 The coupling constant g ~ 0.48eV has been adjusted to yield the best fit. Note 
that in this case To, Eq. (|32|), is of order 20K and hence one observes a very small negative intercept 
of the resistivity, ~ 20/if2cm, as seen in experiment. Inset: the Hall resistivity as a function of 
temperature in comparison with experiment for the same system. 

FIG. 18. The sensitivity of p X x(T) with respect to the band and incommensuration parameters. 
The solid line shows the same p X x(T) as in Fig. |l7| with g = leV. The dashed and dotted lines 
show p X x(T) at a n = 30% doping level and zero incommensuration (AQ = 0) respectively. 
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FIG. 19. The sensitivity of p xx (T) with respect to the changes in spin-fluctuation parameters. 
The curves correspond to (top to bottom) uj sf = 20 + 0.2T(K), 60 + 0.2T(K), 20 + 0.6T(K) and 
60 + 0.6T(K) respectively. All other parameters are assumed the same as in Fig. [l8|. Inset: a xy (T) 
for the input parameters defined in the figure. 

FIG. 20. (a) p xx (T) in the PG regime: the parameters used are: £ _1 = 0.1 + T/1000(K), 
uj s fC = 50meV, g = leV and n = 15%. The solid, dashed and dotted lines correspond to 
—t'/t = 0.2, 0.3 and 0.4 respectively. The effective interaction is assumed to have incommen- 
surate peaks at Q + (±AQ,0) and Q + (0,±AQ), where AQ = 0.25-7T. The three curves show 
a vast difference in the crossover temperature Tq(T). (b) p xx (T) in comparison with experiment 
(Ref. 36). The solid line corresponds to t' = — 0.3i, with g adjusted to yield the same slope as seen 
experimentally. The dashed line corresponds to t' = — 0.25i, with the FL scattering included (see 
text). 

FIG. 21. Comparison of the calculated value of cot &h with experiment (Ref. 36). The solid 
(dashed) line corresponds to the result depicted by the solid (dashed) line in Fig. [20|b, 

FIG. 22. A model FS for a system with strong a-b plane anisotropy. The anisotropy parameter 
is assumed to be a = 0.5 (see Eq. (^)), the doping level is n « 15%. 

FIG. 23. The resistivities along a and b crystallographic directions as a function of temperature 



for the system with the FS depicted in Fig. 22J, using the spin fluctuation parameters appropriate 
to YBa2Cu408 system. Note that there is not only a large quantitative difference, but a qualitative 
difference as well, as is seen in experiment. The anisotropy of the distribution of the hot regions 
leads to the anisotropy of the resistivities. 

FIG. 24. The Hall conductivity as a function of temperature for the system depicted in Fig. 



23. Note that, as usual, a xy reflects an average contribution of the hot and cold regions and hence 



is proportional to T 3 as in the isotropic case (see Fig. 17 and its inset) 
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